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ABSTRACT 

We describe a new radial velocity survey of T Tauri stars and present the 
first results. Our search is motivated by an interest in detecting massive young 
planets, as well as investigating the origin of the brown dwarf desert. As part of 
this survey, we discovered large- amplitude, periodic, radial velocity variations in 
the spectrum of the weak line T Tauri star LkCa 19. Using line bisector analysis 
and a new simulation of the effect of star spots on the photometric and radial 
velocity variability of T Tauri stars, we show that our measured radial velocities 
for LkCal9 are fully consistent with variations caused by the presence of large 
star spots on this rapidly rotating young star. These results illustrate the level 
of activity-induced radial velocity noise associated with at least some very young 
stars. This activity-induced noise will set lower limits on the mass of a companion 
detectable around LkCa 19, and similarly active young stars. 

Subject headings: stars: pre- main sequence — stars: activity — (stars:) plane- 
tary systems: formation — stars: individual (LkCa 19) 
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Introduction 



The first successful detection produced by radial velocity surveys of a planet orbiting a 
star occured in 1 995 when a near Jupite r-mass planet was found around the main sequence 
star 51 Pegasus (IMayor &: Queloall995l ). Subseq uent observations verified the existence o f 
the companion and refined its orbital parameters (iMarcy fc Butlerlll995l ; iMarcy et al.lll997l ). 
By 2007, over 200 planetary systems have been discovered via the radial vel ocity method 



and the statistics of various properties of the planets have been catalogued (IButler et al. 
2006h 



Most radial velocity surveys have been conducted on main sequence stars typically older 
than 1 Gyr. The use of the radial velocity method on younger stars is challenging because 
of increased activity, specifically in the form of signals caused by star spots. Spots on the 
surface of the star rotate in and out of view of the observer, potentially creating spurious 
radial velocity signatures. 



Paulson et al.l (12004 ) looked for planetary companions to stars in the Hyades cluster 



^790 Myr). They found activity induced jitter c onsistent with the empi ric al relationship 



betwe en activity and radial velocity described in ISaar fc Donahue! (119971 ). iPaulson et al. 



( 120041 ) also found one binary [K = 1152 m s x ) and two ca ndidate binaries with li near 



velocity trends, but no planetary companions. More recently, IPaulson fc Yeldal (120061 ) ex- 
amined 61 stars in the (3 Pic moving group (~12 Myr) and the Ursa Majoris association 
(~300 Myr) and found no Jupiter mass planets in short period orbits. The authors found 
these results unsurprising, as planetary companions are rare around main sequence stars 
with the small orbital separations (< 0.1 AU) that these surveys were able to probe. They 
measured maximum radial velocity amplitudes of 30-100 m s _1 for targets in Ursa Majoris 
and 250-600 m s _1 for targets in /3 Pic. 

These searches were motivated, as is ours, by an interest in understanding the formation 
of both stellar and planetary systems. By probing the youngest stars for planetary compan- 
ions, w e can learn at what stage planets form. Since the time scale for the core accret ion 
model (jPollack et al.lll996l ) is longer than formation due to disk instability (jBossl 119971 ). a 
detection of a young planet around a T Tauri star should also provide clues to the formation 
mechanism itself. We can also explore the level of stellar activity in young stars and how 
this can e ffect our ability to me asure radial velocity variations, and therefore detect planets. 
Recently, ISetiawan et al.l (120071 ) have reported the detection of a planet around the 100 Myr 
star HD 70573. 



An additional motivation to search for companions around young stars is to understand 
the origin of the so-called "brown dwarf desert." As more planetary companions were dis- 
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covere d in the last decade, and the completeness of the surveys increased, iMarcy fc Butler 
( I2000l ) noted the dearth of substellar objects with masses greater th an 20 M. Tlip in close or- 
bits (< 3 AU). This "brown dwarf desert" was further quantified by iGrether fc Lineweaver 
(120061 ) who analyzed an unbiased sample of nearby sun-like stars to measure the rate of 
stellar, brown dwarf, and planetary companions. They found that ~16% of their sample 
stars had close (P < 5 yr) companions with masses greater than that of Jupiter. Of these, 
11% ± 3% were stellar, < 1% were brown dwarfs, and 5%± 2% were giant planets, where a 
brown dwarf is defined as an object in the mass range between 13 M Jup to 80 M Jup . 



Armitage fc Bonnelll (120021 ) suggested that the kinematic viscosity of a sufficiently mas- 



sive circumstellar disk could cause brown dwarfs in close orbits, interacting dynamically with 
the disks, to catastrophically migrate towards the host star and be lost to mergers. They 
predict that young stars should have an order of magnitude more brown dwarfs in close 
orbits than main sequence stars, as the time scale for catastrophic migration is ~ 1 Myr. A 
search for substellar companions around T Tauri stars could confirm or refute this scenario, 
improving our understanding of the brown dwarf desert in older stars. 

In this paper, we summarize our spectroscopic and photometric observational techniques 
for our sample of T Tauri stars and report a detection of periodic radial velocity variations in 
the weak line T Tauri star LkC a 19, for which we have a large number of obser vations. LkCa 
19 has a vsini = 18.6 km s -1 (jHartmann et al.l 119871 ) and sp ectral type K0 (IHerbig &: Bell 



19881 ). LkCa 19, with its short rotation period of 2.24 days (IBouvier et al.lll993l ). provides 
an excellent example of the effects of star spots on radial velocity measurements. Using new 
simulations, we show that the radial velocity properties of LkCa 19 are well modeled by the 
presence of one large spot. 



2. Sample, Observations, Data Reduction 
2.1. Sample 

Our total sample consists of 100 young stars in the Taurus star forming region and the 
Pleiades open cluster, with a V magnitude range from 9 to 15. Approximately one third 
of the Taurus objects are classical T Tauri stars, with the remainder weak-line T Tauri 
stars. Stars in the Pleiades and a subsample of our Taurus objects were observed as part 
of the Space Interferometry Mission (SIM) Planetquest precursor science program. These 
observations will identify Jupiter-mass and greater spectroscopic companions; these systems 
will be excluded from the SIM target list. 

Most targets were observed at least once per observing run, weather permitting. A 
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subset of targets were selected for observation on each night of every observing run. Since 
rotation periods for T Tauri stars are typically 10 days or less, this was done to quantify the 
radial velocity variability of a sub-sample, and to understand the possible sources of error 
and velocity jitter such as star spots or other sources of activity (e.g. accretion). 



2.2. Observations : Spectroscopy 



In 2004, we began a radial velocity survey to search for planets a nd brown dwarfs 
around a sample of T Tauri stars using the Coude echelle spectrograph (ITull et al.l Il995l ) 
on the 2.7m Harlan J. Smith telescope at McDonald Observatory. The spectrograph yields 
R = = 60,000 with a 1.2 arcsecond slit. Observations at the 2.7m Harlan J. Smith 
telescope were taken in November of 2004, December of 2004 to January of 2005, November 
of 2005, and February 2006. The Julian dates of these observations are listed in Table 1. 
Current high precision radi al velocity surveys routinely achieve velocity precision of ~ 3 m 



(e.g. 



Butler et al.l 119961 ) . usually by the use of an iodine cell to superimpose calibration 
absorption lines onto the stellar spectrum. However, because we expect T Tauri stars to 
have intrinsic radi al velocity noise considerably higher than th is due to photospheric actvitiy 
such as spots (e.g. IPaulson et al.ll2004J : IPaulson fc Yeldall2006l ). we did not require such high 
precision measurements. We make use of thorium-argon lamps taken before and after each 
object. This method also avoids contaminating the spectra with iodine absorption lines. 
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2.3. Observations : Photometry 



We obtained near-simultaneous differential CCD photometric observations of LkCa 19, 
beginning in December of 2004. The purpose of these observations was to look for any 
concurrent photometric variability and to see if any signal in the photometry is also present 
in the radial velocity variations. Photometric variability in weak line T Tauri stars such as 
LkCa 19 is believed to originate in dark star spots moving in and out of view as the star 
rotates. The resulting false radial velo city modulation sho ws a quarter-phase offset between 
the photometric and RV light curves (jOueloz et al.ll200ll ). A substellar companion around 
a T Tauri star should not produce a photometric signature, assuming there is no planetary 
transit. 



Photometry was obtained using the prime focus corrector (jClaverl Il995t ) on the 0.8 m 
telescope at McDonald observatory. Our observations were taken through Bessel V, B, and R 
photometric filters. The field of view for the images was 42.6 x 42.6 arcminutes. Integration 
times ranged from 3 to 300 seconds. 



2.4. Data Reduction 

2.4-1- Spectroscopy 

All data reduction, including wavelength calibration was performed using standard 
echelle spectral reduction routines in IRAF. The raw spectra were bias subtracted and flat 
fielded using the normalized spectrum of an internal lamp. Optimal extraction, to remove 
cosmic rays and improve signal, was used on all targets except for the bright radial velocity 
standards. The dispersion solution was derived from the time weighted average of thorium- 
argon lamps taken both before and after each image. Individual thorium-argon lamps had 
lines identified visually and calibrated using the NOAO spectral atlas database. We used 7 
to 10 hand-chosen lines per order across ~ 15 orders of the total spectrum for the dispersion 
solution. The final fit was done in both the dispersion and cross-dispersion axis, and has a 
typical RMS value of about 0.002 pixels, or 4 m s _1 , which is a negligible source of error 
for this program. 

The RVSAO package in IRAF was used to determine the radial velocities. IDL was used 
to further analyze and visualize the radial velocity data. The xcsao task in the RVSAO 
package was performed on the extracted and wavelength calibrated spectra to determine 
the radia l velocities. Th i s rou tine is a fourier-based cross-correlation program, described in 



detail by lKurtz fc Mink! (119981 ). Cross-correlation is a powerful tool when coupled with large 
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Table 1. LkCal9 Radial Velocities 



m/s 



JD 




RV 




Error 


S/N 


2453367, 


.8252 


2193. 


.6 


190. 


.5 


81.6 


2453369, 


6507 


1995. 


.0 


141. 


.2 


87.6 


2453370. 


.8191 


540. 


.5 


140 


.1 


52.0 


2453371. 


.6180 


1381. 


.4 


147. 


.9 


76.1 


2453372. 


.6210 


1092. 


.6 


153 


,4 


27.5 


2453373. 


7105 


1095. 


.8 


175. 


,4 


77.6 


2453374. 


7654 


1213 


.0 


331 


.5 


35.4 


2453375. 


.6835 


684. 


.5 


186. 


.8 


59.7 


2453693 


.8006 


1780. 





168. 


.4 


54.2 


2453694. 


.8526 


1193. 


.0 


208. 


.9 


58.4 


2453695. 


.7514 


1160. 


.2 


202. 


.6 


71.3 


2453696. 


.7825 


1555. 


.0 


140 


.5 


81.5 


2453697. 


7609 


451. 


.8 


151 


.4 


78.5 


2453769, 


5967 


309. 


.0 


140. 


.4 


66.5 


2453769. 


7395 


803. 


.8 


158. 


.5 


58.5 


2453770. 


,5880 


2041, 


.2 


195 


.8 


60.8 


2453770. 


7327 


1924. 


.2 


182 


,8 


69.2 


2453771. 


6556 


192 


.6 


280. 


.6 


54.0 


2453771. 


8006 


505. 


.7 


140. 


.4 


46.6 


2453772. 


5921 


2052. 





163 


.2 


57.8 


2453773. 


.5887 


105. 


.4 


144. 


.7 


64.5 


2453773. 


.7526 


0. 


.0 


121 


,0 


57.7 


2453775. 


.6523 


743. 


.8 


148. 


.4 


57.7 


2453775. 


.7981 


307. 


.7 


164. 


.7 


63.2 
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spectral ranges to determine precise radial velocities. We performed cross- correlation using 
several echelle orders, with each order covering approximately 100 angstroms. In order to 
avoid spectral type mismatching, we used the target star itself as the template for the cross 
correlation. Therefore, our velocities are measured relative only to one observation epoch. 
A spectrum with relatively high counts and minimal cosmic ray contamination was chosen 
to be the template. 



2.4-2. Photometry 

Photometric data reduction followed standard procedures using IRAF including bias 
correction, overscan correction, and flat fielding. Stellar objects are identified with the 
DA0PH0T package within IRAF and then aperture photometry was used on all stars in the 
field. The error in the photometry was estimated by selecting a background star, or several, 
and looking at the standard deviation of that star over all observations. Such stars should 
show no photometric variability, and any scatter in their measurements provides an easy 
estimate on the accuracy of the photometry. This error is ~0.01 magnitude. 



3. Analysis 
3.1. Radial Velocity Uncertainties 



In order to gain a better understanding of the accuracy of our rad ial velocity technique , 
we observed s ix rad ial ve locity standard s t ars, s electing objects from lNidever et al.l (120021 ). 
Butler et al.l (119961 ). and ICumming et al.l (119991 ) . By measuring radial velocities of these 
targets each night of each observing run, we estimate the systematic errors and the accuracy 
of the radial velocity measurements with our instrumental setup. We observed three bright 
standard stars, in the 3rd to 4th V magnitude range, in addition to three fainter objects 
around 6th to 7th V magnitude. These standard stars were reduced and analyzed in the same 
way as our target stars, using one epoch of the standard star itself as the cross-correlation 
template. Within a given run, we obtained very low radial velocity scatter, however there 
are systematic effects between runs that appear in the observations of our standard stars, 
[t!] 

Figure [T] shows RV measurements for our 3 bright standard stars. Each standard shows 
a distinct systematic offset between the several observing runs. Simple explanations do not 
account for this systematic error - the residuals in the dispersion solution look very similar 
in each run and we accounted for the motion of the earth itself. We suspect the likely cause 
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Fig. 1. — Radial velocity shifts of the standards 107 Psc, Tau Ceti, and HD 4628 showing 
systematic offsets between observing runs. 



of the offset is the resetting of the slit plug from one observing run to the next, which could 
alter the position of the slit relative to the optical axis of the telescope and the thorium- 
argon lamp. However, even with this systematic error uncorrected, the overall scatter of the 
standard stars remains ~ 120 m s _1 (see Table 2), which we believe is our nominal level of 
precision. The final error listed in Table 1 is the sum in quadrature of this systematic error 
and error estimated from the order-to-order scatter in the measured radial velocities of LkCa 
19. 



3.2. Detection of Periodic Signals 

Our T Tauri targets have considerably larger radial velocity variations, with standard 
deviations in the range of ~ 200 - 3000 m s _1 , than the radial velocity standards. LkCa 19 
shows radial velocity variation with a standard deviation of ~ 670 m s -1 , which is about six 
times higher than our nominal radial velocity standard. 

For targets with a significant number of observations, such as LkCa 19, we analyzed the 
data to look for any periodic variations of the radial velocity measurements. This w as done 



with an IDL implementation of the Scargle method of power spectrum estimation (jScargle 



19821 ). This method is best suited for our data, as it does not require the data to be evenly 



spaced in time. The LkCal9 power spectrum is shown in Figure [2j 
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Table 2. Radial Velocity Standards Observed 



Star 


N 


V Mag 


Standard Deviation (m/s) a 


107 Psc 


20 


5.2 


100 


HD 4628 


20 


5.75 


129 


Tau Ceti 


19 


3.50 


184 


HD 65277 


15 


8.12 


110 


HD 80367 


19 


8.16 


107 


HD 88371 


18 


8.43 


132 



Includes systematic error seen in Figure [TJ 
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1 2 3 4 5 6 



Period (days) 

Fig. 2. — The Scargle method power spectrum on all radial velocity points for LkCal9. The 
peak with the largest power corresponds to a false alarm probability < 0.001 

Once the initial power spectrum was computed over a wide range of periods, the power 
spectrum was computed again with dense sampling in order to resolve the strongest peak 
and to determine the best period, which we found to be 2.177 days with an uncertainty of 
0.0264 days. We then used the period to phase fold the original data so that any radial 
velocity periodicity could be evaluated (Figure [3]). 

Further testing of the power spectrum analysis is necessary to ascertain the true signifi- 
cance of the detection of a period. This was accomplished with a Monte Carlo simulation of 
our radial velocity data. For LkCal9, 100,000 simulated sets of normally distributed random 
radial velocities with identical overall velocity scatter to our observed data were analyzed. 
These random data sets were analyzed using the same temporal cadence as our actual ob- 



- 11 - 



2500 




T3 

(2 500 



Phase 

Fig. 3. — Phase- folded LkCa 19 radial velocities with a period of 2.177 days. The period 
was determined from a Scargle method power spectrum. Error bars are a combination of 
the standard deviation of the mean of the six echelle orders used to compute radial velocity 
added in quadrature with the systematic errors from the radial velocity standards of 120 
m/s. No correction for the systematic errors between runs was applied. 



servations. None of these instances of random velocity data had a peak as strong or stronger 
in its power spectrum analysis than is present in the actual LkCa 19 data. This indicates a 
high significance of this period detection, and a false alarm probability (FAP) of less than 
10~ 5 . Thus, we conclude that LkCal9 shows periodic radial velocity variations with a pe- 
riod of 2. 177±0.0264 days. H owever, this period is very near the 2.24 day rotation period of 
LkCa 19 (IBouvier et al.lll993l ). suggesting that the observed signal is not the result of orbital 
motion. Additional periods in the power function occur near 2.238 and 2.251 days, which 
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are near the power of the strongest peak at 2.177 days. This is an additional indication 
that we our observing the rotation period of the star. The application of a correction for the 
systematic effects seen in Figure 1 produced no noticeable change in the Scargle periodogram 
of LkCal9, and no change in any of the detected periods. 



3.3. Bisector Analysis 



Bisector analysis of line profiles was first suggested by lGrayl (119761 ) as a means to quantify 
the asymmetry of a specific absorption or emission line. Gray's interest was in the intrinsic 
line profile shapes of stars of various spectral types, and as a probe of granulation due to 
convection. Bisector analysis as a means to check on the origin of radial velocity vari ations 
started with analysis of the original extrasolar pl a net 5 1 Pegasus (IHatzes et al.lll997l ). The 
method was used subsequently by iQueloz et al.l (120011 ) to look for a correlation between 
the slope of the bisectors of the spectrum's cross corre lation function and the measured 
radial velocity of the star. Similar analysis was done by iMartmez Fiorenzano et al.l (120051 ) 



on several potential extrasolar planets, and by iBouvier et al. J2007h to test the origin of 
apparent radial velocity variations in the T Tauri star AA Tau. We performed bisector 
analysis on the cross correlation function (CCF) used to measure the radial velocity shift 
from one observation to another to test whether spots might be the cause of the observed 
radial velocity variations in LkCa 19. 

To verify that we do not see correlated radial velocity and CCF bisector changes in 
systems with companions, we ana lyzed our data on a known extra solar planet host from 
the Keck High Precision Survey (jVogt et al.ll2002l ). HD 68988 was observed to test our 
ability to measure a low amplitude periodic signal. The "span" of the bisector is similar to 
the slope of a line; it is measured by comparing the values of the bisector at two distinct 
heights of the cross correlation function (Figure H|). Bisector analysis of HD 68988 yielded 
no correlation between the bisector span and radial velocity, consistent with its classification 
as a planet-hosting star. However, analysis of the bisectors for LkCa 19 through the most 
recent observing runs indicate a strong correlation between the bisector span and the radial 
velocity, a clear indication of star spot-induced radial velocity variations. These results are 
shown in Figure [51 

It is possible to see the variation in line profile shape in the LkCa 19 spectra. In order to 
maximize signal to noise, we co-added several nights of similar velocity from our observations 
in February of 2006. In addition, the averaged spectra were continuum normalized and 
smoothed with a moving average of 3 pixels. The shifting asymmetry of these line averages 
can be seen in Figure 
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Radial Velocity (km/s) Radial Velocity (km/s) 

Fig. 4. — Top Left: A cross-correlation function and its bisector. The dotted line connects 
opposite sides of the CCF curve - the midpoint of this line represents one of the points of 
the bisector. Bottom Left: Two mean velocities are computed - one in the upper region 
denoted by the dashed lines and another between the dotted lines. The bisector span is the 
difference between these two velocities. Right: Four bisectors for LkCal9 at opposite phases 
of the radial velocity curve. The changing asymmetry of the profile is noticeable. Dotted 
lines indicate the CCF power ranges in which the span is measured. 



3.4. Photometric Analysis 

Photometric observations were obtained nearly simultaneously with our spectra of LkCa 
19 to further aid the interpretation of any radial velocity variations. The photometric light 
curves and any periodicity can be compared to the radial velocity curves and detected peri- 
ods. 

A FORTRAN code provided by J. Anderson was used to analyze the multiple exposures 
over many nights and to create light curves for our target stars. The photometric zero-point 
was determined using a sigma-clipped average of all the stars in the field. The photometry 
was phase-folded using the radial velocity period found in §3.2. Furthermore, the photometry 



-14- 




Fig. 5. — Left: LkCal9 bisector analysis. The radial velocity and bisector appear to be 
correlated, with a linear correlation coeffi cient of —0.7 3. The slope is consistent with the 
spot induced radial velocity variations of iQueloz et al.l (120011 ). Right: HD 68988 bisector 
analysis. No correlation is found. 



was analyzed to look for other periodicity. 

We analyzed 41 R band photometric points for LkCa 19, which are binned into nightly 
averages producing 17 photometric measurements. On any given night, the R-band images 
were taken within ~20 minutes of each other. The data points are plotted, phased to the 
radial velocity period of LkCa 19, in Figure [71 The photometry does not appear in phase 
with the radial velocity signal. A Scargle power spectra analysis for our photometry yields a 
period of 2.05575 days. However, this photometric period determination does not withstand 
Monte Carlo test scrutiny, and the false alarm probability for the 2.05575 day LkCa 19 
photometric period is 40%. The poor strength of this period is likely due to the smaller 
number of measurements compared to the radial velocity data, and the greater noise in the 
photometric measurements compared to the amplitude of the photometric variability. This 
period, however, is close t o both the detected radial velocity period and the known 2.24d 



rotation period of LkCa 19 (IBouvier et al.lll993l ). The peak-to-peak amplitude in the R-band 



is 0.12 magnitudes. We observed the same amplitude in V-band, however in the B-band the 
peak-to-peak amplitude is larger at 0.18 magnitudes. 
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Fig. 6.— LkCa 19 line profiles of the He I line at 6678 A. The "high" and "low" radial 
velocity averages were created by averaging observations of similar radial velocity, 3 spectra 
from the high velocity case and 4 from the low velocity case, from consecutive nights of 
observing run in February of 2006. The spectra were heliocentric corrected before averaging 
and smoothed with a moving average of 3 pixels after they were combined. 



4. Star Spot Models 

4.1. Conceptual Framework 

The goal of our star spot modeling was to test the hypothesis that at least some of the 
radial velocity signatures observed in T Tauri stars can be attributed to the star's rotation 
combined with the presence of low temperature spots on the surface of the star. In order to 
test these assumptions, it is necessary to numerically replicate the effect of star spots on a 




Fig. 7. — LkCal9 RV and R-band photometry phased to periods of 2.177 and 2.05575 days. 



synthesized stellar spectrum. 

To accomplish this task, a stellar disk integration model was constructed. The apparent 
stellar disk is divided into sections, allowing for portions to be at a lower temperature 
and therefore lower intensity than the rest of the star, with appropriately different spectral 
features. With the modeled stellar disk divided between a stellar spectrum at one effective 
temperature and one or more localized spots with a lower temperature, an integrated spotted 
star spectrum can be created and analyzed in the same manner as the observed stellar data. 
This modeling involved integrating an artificial stellar disk, with control over the location 
and size of spots on the apparent disk, and was implemented in the IDL language. 



In order to construct the apparent stellar disk of a star containing spots, it was necessary 
to first synthesize appropriate spectra for both the star and t he spot. The input spectra were 
produced using the SME (Spectroscopy Made Easy) code of IValenti fc Piskunov! (119961 ). We 
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chose a 100 A bandpass near 6300 A, similar to one of the echell e orders in our actual data . 



The SME synthesis code utilized NextGen atmospheric models (lAUard &: Hauschildtlll995l ) 
and a line list from the Vienna Atomic Line Data (VALD) database. The dispersion of the 
artificial spectra was 0.02 A per pixel, with the input and output spectra having 5000 pixels. 
The stellar temperature was chosen to be 4300K and the spo t temperature 2900K, based 



on observations and modeling of V410 Tau (IPetrov et al.lll994l ). This produced a brightness 



ratio between the spots and the stellar surface of 0.08. 

Spots were defined by the colatitude and longitude of the spot center, and the radius of 
the spot. The spot's projected shape and location on the stellar disk were used to trigger the 
integration routine to use the lower temperature and intensity spot spectrum rather than the 
standard stellar spectrum when the given portion of the stellar disk is included. Each grid 
point of the star that was deemed to be inside the spot was assigned the cooler, dimmer, 
spot spectrum. Additionally, each grid point has the corresponding radial velocity shifts 
based on the rotation of the star. Therefore, Doppler broadening of the line profiles was 
handled within the disk integration itself. Stellar limb darkening was handled by using the 
appropriate SME spectra for multiple limb angles. The disk integration's output was tested 
against the disk integration output of the SME code - a star with no spots - and correctly 
returned a spectrum identical to the output of the original SME disk integration. 



4.2. Analysis of Simulated Spectra 

The simulation produced a series of spectra, which were then read into IRAF and 
analyzed for radial velocities in a manner identical to the actual observations. Because the 
relative spectral fluxes are preserved, the IRAF bandpass spectrophotometry task sbands 
was used on a continuum region of the spectrum to estimate the photometric variations of 
the star due to the spots, using a wavelength bandpass of 1 A wide centered at 6235.5 A. 

Twenty synthesized spectra were created representing 20 evenly spaced rotation steps 
of 18 degrees and thus one complete stellar rotation. Each spectrum was cross-correlated 
against a template synthetic spectrum of the same vsiiai but without any star spots. The 
primary difference in the model analysis was that, without multiple echelle orders, we could 
not estimate the radial velocity error from the order to order scatter. The artificial spectra 
were otherwise run through identical analysis as the actual data, although they were not 
rebinned or resampled to the resolution of the observations, nor was any noise added. The 
RVSAO package was used to compute the velocities of the simulated spectra, and the bisector 
analysis was performed in the same manner as for the actual stellar data. 
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The simulation code can be used with a variety of stellar spot configurations, ranging 
from a simple equatorial spot in an uninclined star, to a polar spot on a star with a tilt, to 
a multispot configuration, or any combination of the above. This wide range of possibilities 
led to the question of which scenarios to simulate, and if any of them could reproduce the 
magnitude of radial velocity variations seen in our program stars. We discuss three specific 
cases. All of the computations were done with a vsini value of 20 km/s, which is near the 
value for LkCal9, unless otherwise noted. 



4-2.1. A Circumpolar Spot 

A single equatorial spot would yield a flat radial velocity curve while the spot is on 
the opposite side of the star from the observer. In order to produce, with a single spot, a 
radial velocity curve that shows persistent amplitude modulation over time, it is necessary 
for the spot to always be visible to the observer. Therefore, a spot close to the pole of 
the star, with the star inclined towards the observer, has been chosen as a simple geometry 
to model. In this case, the spot has a colatitude of 40° and a radius of 15°; the stellar 





Strassmeier et al. 


1994; 


J one our et al. 


1994; 


Johns-Krull & Hatzes 


1997: 


Unruh et al. 


2004 


). The circumpolar spot rotation sequence is illustrated in Figure [U 



The photometric light curve and the radial velocity curve are shown in Figure [H This 
spot configuration produces a radial velocity curve that varies continuously and has an 
equally regular photometric curve. For this spot the photometric variation is at most 0.086 
magnitudes and the full radial velocity a mplitude is 1820 m/s . We adopt a maximum flux- 
weighted filling factor for the spot following lO'Neal et all (119961 ). who define the flux weighted 
filling factor as the total fractional projected area of spots on the observed hemisphere 
weighted by limb drakening. This is readily comput ed via our disk i n tegrat ion code and 
is 9.1% comp ared to the predicted values of 8.7% of ISaar &; Donahue! ( 119971 ) and 6.4% of 
Hatzed (120021 ) . Both papers contain an empirical relationship between the related filling 
factor, radial velocity amplitude, and vsini. 

Bisector analysis of this spot configuration (Figure [9]) shows the same general trend as 
the bisector analysis of LkCa 19. The width of the linear bisector span - velocity correlation 
appears intrinsic, not merely an artifact of noise in the actual stellar data. The correlation 
between the span and the measured velocity is degenerate, in that the same velocity does 
not always produce an identical bisector span, even for one unique spot configuration. 



-19- 




Fig. 8. — The rotation sequence of a single spot rotating around the pole of an inclined star; 
i = 35° and the spot's colatitude is 40°. The X marks the pole. 



4-2.2. Fifteen Random Spots 

In addition to a simple one spot model, we performed a simulation that involves many 
spots, randomly distributed, with the star at a random inclination. The random spot distri- 
bution was accomplished using the randomu task in IDL to produce 15 uniformly randomly 
distributed colatitude and longitude pairs. The sizes of the spots were fixed with five spots 
each of three different sizes (Figure [TO]) . 
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Fig. 9. — Radial Velocity and Photometric curves for the single circumpolar spot (left), and 
the corresponding bisector span vs. radial velocity plot (right). 




Models with multiple spots of different sizes produced very different radial velocity and 
photometric light curves (Figure [TT1) than the simpler one spot model. The bisector analysis 
demonstrates the intrinsic scatter of the bisector span - radial velocity correlation. There 
is a degeneracy in the bisector correlation such that the bisector alone can not predict the 
measured velocity. The trend observed in this simulation is consistent with the correlation 
seen in LkCa 19. The most significant feature of this simulation is the amplitude of the radial 
velocity curve, roughly 1200 m/s, and the amplitude of the photometric measurements, ~ 
0.1 magnitudes. 
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Fig. 11. — Radial velocity and photometric curves for fifteen random spots (left), and the 
corresponding bisector span vs. radial velocity (right). The correlation coefficient of a linear 
regression fit is -0.91, indicating a strong correlation. 



4.3. The Relationship between Photometric and Radial Velocity Amplitude 

The results from the three example spot simulations, along with other spot geometries 
not described here in detail, prompted an inquiry into the relationship between the pho- 
tometric amplitude and the radial velocity amplitude. We compared the photometric and 
radial velocity amplitudes of several simulations, along with our actual observations. In- 
cluded were four circumpolar spots of various sizes similar to the spot sequence described in 
§4.2.11 Similarly, several additional spot simulations that contained 15 random spots were 
run, each with a different set of spot locations, but with the same inclination and spot size 
distribution. The results for the simulations and LkCal9 are plotted in Figure [T2] In all 
cases v sinz = 20 km s _1 . 

LkCa 19 falls very close to the line representing circumpolar spots of various sizes. 
The multi-spot configuration has a distinctly different slope. The LkCa 19 simulation point 
in Figure [15] represents a best guess at the parameters of LkCa 19 itself. The inclin ation 



angle was estimated at 25° based on the measured vsini value (IHerbig fc Belli I1988T ) and 
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Fig. 12. — Full Radial Velocity amplitude versus full photometric amplitude for a variety of 
simulations, observations of LkCa 19, and our simulation of LkCa 19. 

an inferred radius based on its luminosity and temperature. The actual vsmi value of 18.6 
km/s was also used in this simulation, otherwise the parameters are similar the circumpolar 
spot model described earlier. 



Table 3. Correlation Statistics for Bisector Span and Radial Velocity 



Object 


Spearman p a 


Two Sided Significance 


Linear R 


LkCal9 


-0.79 


3 x 10~ 5 


-0.73 


Simulation (15 Spots) 


-0.94 


1.4 x 10~ 9 


-0.93 


Simulation (One Circumpolar Spot) 


-0.96 


8.6 x 10~ 12 


-0.92 



a Rank Statistic 



-24- 



5. Conclusion 

Analysis of early results from our radial velocity survey of young stars reveals large 
radial velocity variations in our target stars compared to radial velocity standards. LkCa 
19 shows a radial velocity signal with a significant periodic variation. Initial analysis of the 
radial velocities of LkCa 19 led to twice-nightly observations in our February 2006 run to 
obtain better phase coverage. Bisector analysis, however, indicates that the radial velocity 
variations in LkCa 19 are very likely induced by star spots, as does the similarity between 
LkCal9's 2.177±0.02 64 day radial velocity period and its known rotation period of 2.24 days 



flBouvier et alJll993 ). 



While it is possible that a close planet could experience tidal locking (IMarcy et al.lll997l ) 
and thus have an orbital period identical to the rotation period of the star, it is more likely we 
are seeing radial velocity modulation arising in star spots. A surprising feature of LkCa 19 is 
the phase coherence of the spot-induced radial velocity modulation. The velocities remain in 
phase at the period of 2.177 days over a time span of approximat e ly 13 months. This result 



is similar to the false planet detection reported by iQueloz et al.l (120011 ). also attributed to 



star spots, which showed reasonable phase- coherence for a time span of roughly 2 years. 

We have also presented a new simulation of the effects of a rotating star with spots 
on line profiles and have measured the radial velocities and photometric variations from 
spectra produced by the simulation. The models can reproduce the photometric and velocity 
amplitudes of LkCa 19. All simulations show a strong and significant bisector correlation, 
summarized in Table 3. Further, each spot simulation shows a 0.25 lag in phase between 
photometry and radial velocity, regardless of spot configuration, which is expected for cool 
spots. 

The spot geometry with the most regular or planet-like radial velocity curve is that of 
a single spot near the pole of a star with a low inclination angle. A highly spotted star with 
the spots set randomly on the stellar surface, produces much more irregular radial velocity 
and photometric light curves. 

If the magnitude of the photospheric activity induced radial velocity variations of LkCa 
19 is typical of young stars, then our ability to find companions to these objects will be 
limited to the extremely massive and short-period sub-stellar companions. LkCal9 has an 
exceptionally short rotation period and relatively high v sin i value, which offers some hope 
that its spot-induced radial velocity variations are not necessarily typical for T Tauri stars. 



We made use of IDL libraries available from the NASA Goddard Space Flight Center, 
and the John Hopkins / Advanced Physics Lab / Space Oceanography group, as well as the 
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IRAF, originally from the National Optical Astronomy Observatory, and now maintained by 
the volunteers at IRAF.net. Jay Anderson provided very helpful assistance with the pho- 
tometric reduction, and useful discussions regarding the spot simulation code. The authors 
are grateful to the anonymous referee for useful comments. This research also made use 
of NASA's Astrophysics Data System. M.H. was partially supported by NSF Cooperative 
Agreement Number HRD-0450363 as a part of the Rice University Alliance for Graduate Ed- 
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NASA Grant 05-SSO05-0086 to Lowell Observatory (PI - L. Prato), and from the Space In- 
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